The Einstein relation generalized to non-equilibrium 
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The Einstein relation connecting the diffusion constant and the mobility is violated beyond the 
linear response regime. For a colloidal particle driven along a periodic potential imposed by laser 
traps, we test the recent theoretical generalization of the Einstein relation to the non-equilibrium 
regime which involves an integral over measurable velocity correlation functions. 
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A comprehensive theory of systems driven out of equi- 
librium is still lacking quite in contrast to the univer- 
sal description of equilibrium systems by the Gibbs- 
Boltzmann distribution. Linear response theory provides 
exact relations valid, however, only for small deviations 
from equilibrium [![. The arguably most famous linear 
response relation is the Einstein relation 



D = k B Tfi, 



(1) 



involving the diffusion constant D, the mobility //, and 
the thermal energy k B T [2]. In his original derivation 
for a suspension in a force field, Einstein balances the 
diffusive current with a linear drift. The Einstein re- 
lation embodies a deep connection between fluctuations 
causing diffusion and dissipation responsible for friction 
expressed by a finite mobility. 

In the present Letter, we report on the extension of 
the classical Einstein relation beyond the linear response 
regime using a driven colloidal particle as a paradigmatic 
system. Our previous theoretical work [3[ and its present 
experimental test thus introduce a third type of exact re- 
lation valid for and relevant to small driven systems cou- 
pled to a heat bath of constant temperature T. The pre- 
viously discovered exact relations comprise, first, the fluc- 
tuation theorem 0, [1| which quantifies the steady state 
probability of observing trajectories of negative entropy 
production. Second, the Jarzynski relation 0] expresses 
the free energy difference between different equilibrium 
states by a nonlinear average of the work spent in driv- 
ing such a transition Q. Both the fluctuation theorem 
and the Jarzynski relation as well as their theoretical 
extensions || @, have been tested in various experi- 
mental systems such as micro-mechanically manipulated 
biomolecules ffl) . Il2j |. colloids in time-dependent laser 
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traps 
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Rayleigh-Benard convection [lf|, me- 



chanical oscillators [l7(, and optically driven single two- 
level systems [l8j]. Such exact relations (and the study of 
their limitations) are fundamentally important since they 
provide the first elements of a future more comprehensive 
theory of non-equilibrium systems. 

For a non-equilibrium extension of the Einstein rela- 
tion (fTJ), consider the overdamped motion x(t) of a par- 
ticle moving along a periodic one-dimensional potential 



V{x) governed by the Langevin equation 



1 



x(t) = ~F(x(t))+£(t) 

1 



(2) 



with F = —dV/dx + f and / a non-conservative force. 
The friction coefficient 7 determines the correlations 
{£(*)£(*')) = 2(k B T/j)S(t - t') of the white noise f. 
Therefore Eq. ((2]) describes a colloidal bead driven to 
non-equilibrium under the assumption that the fluctuat- 
ing forces arising from the heat bath are not affected by 
the driving. 

For the crucial quantities D and fi, it is convenient to 
adapt definitions which can be used both in equilibrium 
and beyond linear response, i.e., in a non-equilibrium 
steady state characterized by / = const. ^ 0. The diffu- 
sion coefficient is given by 



D= lim [(x*(t)) - (x(t)f]/(2t), 

t—>oo 



(3) 



where (■ • •) denotes the ensemble average. Both theoret- 
ical work [l9( | and a recent experiment (20j have shown 
that the force-dependent diffusion constant can be sub- 
stantially larger than its equilibrium value. The mobility 



/' 



d(x) 



(4) 



quantifies the response of the mean velocity (x) to a small 
change of the external force /. If the response is taken 
at / = 0, which corresponds to equilibrium, one has the 
linear response relation (fT]). How does the Einstein rela- 
tion change for / ^ 0, i.e., what is the relation between 
a force-dependent diffusion constant D(f) and a force- 
dependent mobility Is there a simple relation at 
all? We have recently shown that under non-equilibrium 
conditions the Einstein relation (fT]) has to be replaced 

by 1 



D = k B T^i + dr J(r) 



(5) 



where the second term on the right hand side is given 
by an integral over a known "violation function" I(r) 
involving measurable velocity correlations to be discussed 
in detail below. Such a relation is complementary to 
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FIG. 2: (a) Reconstructed potential V(x). (b) Tilted poten- 
tial. The colloidal particle is subjected to a constant driving 
force / ~ 0.06 pN and the periodic potential V(x). 
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FIG. 1: (a) Experimental setup, (b) Typical trajectory of 
the angular particle position for a mean particle revolution 
time ~ 5.8 s. 



introducing an effective temperature which replaces T in 



Eq. ([T]) in an attempt to keep its simple form [2l|, |2 



It has the advantage that knowledge of I(t) offers us a 
better understanding of the crucial characteristics of the 
non-equilibrium steady state that causes the breakdown 
of the Einstein relation |T]) . 

In our experiment we subject a single colloidal sil- 
ica bead with 1.85 /zm diameter to a non-equilibrium 
steady state by forcing it along a toroidal trap (R = 
1.75 /mi) created by tightly focussed rotating optical 
tweezers [23l Q (see Fig. 1). This is achieved by fo- 
cusing the beam of a Nd:YAG laser (A = 532 nm) with a 
microscope objective (lOOx, NA=1.3) into a sample cell 
containing a highly diluted aqueous suspension of silica 
particles with 1.85 /im diameter. A pair of galvanomet- 
ric driven mirrors with telescope optics deflects the beam 
along a circular path and thus confines the silica bead to 
an effectively one-dimensional motion. Depending on the 
velocity of the rotating trap three different regimes can be 
distinguished (23[. (i) For small velocities friction forces 
are much smaller than the trapping force, the trapped 
particle is able to follow the trap, (ii) With increasing 
velocity the trap is not strong enough to compensate the 
viscous force of the fluid, the particle escapes from the 
laser trap. However, every time the laser passes the parti- 
cle it is still dragged a small distance along the circle and 
moves with a constant mean velocity around the torus, 
(iii) As the focus speed increases (quasi)-equilibrium con- 
ditions are established and the particle is able to diffuse 



freely along the torus. With the trap rotation frequency 
set to 567 Hz the experiments are performed in the in- 
termediate regime (ii) where the particle is observed to 
circulate with a constant mean velocity. Since the dis- 
placement of the particle by a single kick depends on the 
laser intensity and is approximately 10 nm, under our 
experimental conditions the spatial (50 nm) and tem- 
poral (80 ms) resolution of digital video microscopy is 
not sufficient to resolve single "kicking" events. There- 
fore the particle can be considered to be subjected to 
a constant force / along the angular direction x. Ad- 
ditionally the scanning motion is synchronized with an 
electro-optical modulator (EOM) which allows the peri- 
odic variation of the laser intensity along the toroid. In 
the experiment the tweezers intensity P is weakly mod- 
ulated (AP/P < 10%). This small intensity modulation 
superimposes an additional periodic potential V(x) act- 
ing on the particle when moving along the torus. As the 
result, the particle moves in a tilted periodic potential. 
Both the potential V(x) and the driving force / are not 
known from the input values to the EOM but must be 
reconstructed as described in detail below. 

The central quantitity of Eq. ([5|) is the violation func- 
tion I(t) which can be written as Q 



I(r) = ([x(t + r)-(x)}[vs(x(t))-(x)}). 



(6) 



It correlates the actual velocity x(t) with the local mean 
velocity v s (x) subtracting from both the global mean ve- 
locity (x) = 2tt Rj s that is given by the net particle flux 
j s through the torus. In one dimension for a steady state, 
the current must be the same everywhere and hence j s 
is a constant. The offset t is arbitrary because of time- 
translational invariance in a steady state and in the fol- 
lowing we set t — 0. The local mean velocity v s (x) is 
the average of the stochastic velocity x over the subset of 
trajectories passing through x. An equivalent expression 
is j s = v s (x)p s (x) connecting the current with the proba- 
bility density p s (x). The local mean velocity can thus be 
regarded as a measure of the local violation of detailed 
balance. Since in equilibrium detailed balance holds and 
therefore v s (x) = (x) = 0, the violation (J6j) vanishes and 
Eq. (© reduces to Eq. |1|). 

For an experimental test of the non-equilibrium Ein- 
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FIG. 3: a) Experimentally measured violation function I(t) (solid line), b) Comparison of the velocities involved in the 
violation function 7(r). For an ideal cosine potential, we sketch the probability distribution p s (x) (solid gray line), the local 
mean velocity v a (x) together with the drift velocity and their mean (x) versus the angular particle position. The drift velocity 
is the deterministic part F/-y of the actual velocity x. The sign change in 7(r) at (2), (3), and (4) can be understood as follows. 
In a steady state, a single particle trajectory will start with highest probability in the shaded region and, for an illustration, 
we choose its maximum as starting point (1) determining the value v s {x(t)) in Eq. (JSJ). Neglecting thermal fluctuations, the 
particle would follow the dashed line and during a small time step r the product F(x(t + r))v s (x(t)) is positive. If the particle 
passes (2), the product would become negative. The sign changes again if the particle passes (3) and then (4) and so on due to 
the periodic nature of the potential. Thermal noise and averaging over all trajectories does not change this behavior responsible 
for the oscillations of 



stein relation ([5]), we measure trajectories of a single col- 
loidal particle for different driving forces / by adjusting 
the intensity transmitted through the EOM. From a lin- 
ear fit to the data we first determine the mean global ve- 
locity (x). Next, we extract the mean local velocity v s (x) 
from the histogram p s (x) with the coordinate x confined 
to < x < 2nR. Since measurements are performed with 
a sampling rate of 80 ms, we cannot directly access the 
velocity x(t) experimentally. To calculate the violation 
integral I(t), we decompose x(t) into a randomly fluctu- 
ating Brownian part and a drift term, see Eq. @. We 
then transform I(t) as 

I(r) =([v s (x(r)) - (x(t))H(*(0))} - (i) 2 

7 ox (7) 

+ (£(rHOz(0))>. 

The generalized potential &(x) is determined via the 
measured stationary probability distribution, p s (x) = 
exp[— &{x)} [3J. For r > 0, the last term vanishes because 
then £(t+r) and x(t) are uncorrelated. Thus the function 
I{t) depends on two measurable quantities, the current 
j s and the stationary probability distribution p s (x). 

The potential V and the driving force / are determined 
by integrating the force 

dV <9$ 
along the torus. We obtain 



and 

V(x) = k B T<S>(x) + f dx' [f - jv s (x')] (10) 
Jo 

up to an irrelevant constant. In Eq. terms involving 
V and «3> are zero due to the periodicity of our system. 
Both, the potential V(x) and the tilted potential V(x) — 
fx are shown in Fig. 2. The mobility fi = A(x) / Af is 
determined from the change of the global mean velocity 
A(x) upon a small variation of the force Af. 

With the experimentally determined quantities, we 
measure the violation function I(t) shown as solid line in 
Fig. 3a for / = 0.06 pN. It clearly displays the two time 
scales present in the system. First, the driving leads to 
an oscillatory behavior with a period equal to the mean 
revolution time ~ 5.8 s. Second, the diffusion causes a 
broadening of the particle's position resulting in a decor- 
relation between actual and local velocity and hence an 
exponential decay with time constant ~ 2.3 s indicated 
by the dashed line (Fig. 3a). To understand the behavior 
of I(t) in more detail it is helpful to compare the dif- 
ferent velocities involved in the violation function I(t) 
which are sketched in Fig. 3b. 

After numerical integration of the experimentally de- 
termined I{t) we finally calculate the diffusion coefficient 
according to Eq. ([5]). To quantify the relative impor- 
tance of the violation integral we plot the two terms of 
the right hand side of Eq. ([5]) separately for five different 
values of the driving force in Fig. 4. Their sum is in good 
agreement with the independently measured diffusion co- 
efficient directly obtained from the particles trajectory 
using Eq. ([3]). As the maximal error for the indepen- 
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FIG. 4: Experimental test of Eq. ((5]l for different driving 
forces /. The open bars show the measured diffusion co- 
efficients D. The stacked bars are mobility (gray bar) and 
integrated violation (hatched bar), respectively. 



dent measurements we estimated from our data ±3% for 
the diffusion coefficient D, up to ±10% for the violation 
integral, and ±7% for the mobility [i. 

We emphasize that under our experimental parameters 
the violation term dominates the diffusion coefficient (up 
to 80%) and must not be ignored. In Fig. 4 one observes 
a non-monotonic dependence of the violation integral on 
the driving force. This is due to the fact that the maxima 
of n(f) and D(f) do not occur at the same driving force 
but are slightly offset [lfj] ■ This implies for the violation 
function a maximum followed by a minimum as a func- 
tion of /. For very small driving forces, the bead is close 
to equilibrium and its motion can be described using lin- 
ear response theory. As a result, the violation integral 
is negligible. Experimentally, this regime is difficult to 
access since D and /i become exponentially small and 
cannot be measured at reasonable time scales for small 
forces and potentials as deep as 40/cbT (cf. Fig. 2a). 
For much larger forces, the relative magnitude of the vi- 
olation term becomes smaller as well. In this limit, the 
imposed potential becomes irrelevant and the spatial de- 
pendence of the local mean velocity, which is the source 
of the violation term, vanishes. The fact that in our 
regime the violation term is of the same order of magni- 
tude as the mobility proves that we are indeed probing 
the regime beyond linear response. Still, the description 
of the colloidal motion by a Markovian (memory-less) 
Brownian motion with drift as implicit in our analysis 
remains obviously a faithful representation since the the- 
oretical results are derived from such a framework. 



The Einstein relation generalized to non-equilibrium as 
presented and tested here for the driven motion along a 
single coordinate could be considered as a paradigm. Ex- 
tending such an approach to interacting particles and re- 
solving frequency dependent versions of Eq. © [3(] while 
certainly experimentally challenging will provide further 
insight into crucial elements of a future systematic theory 
of non-equilibrium systems. 
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